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The statefinder diagnostic pair is adopted to differentiate viscous cosmology models and it is found 
that the trajectories of these viscous cosmology models on the statefinder pair s — r plane are quite 
different from those of the corresponding non-viscous cases. Particularly for the quiessence model, 
the singular properties of state parameter w = — 1 are obviously demonstrated on the statefinder 
diagnostic pair planes. We then discuss the entropy of the viscous / dissipative cosmology system 
which may be more practical to describe the present cosmic observations as the perfect fluid is just 
a global approximation to the complicated cosmic media in current universe evolution. When the 
bulk viscosity takes the form of f = £id/a(£i i s constant), the relationship between the entropy 
5* and the redshift z is explicitly given out. We find that the entropy of the viscous cosmology is 
always increasing and consistent with the thermodynamics arrow of time for the universe evolution. 
With the parameter constraints from fitting to the 157 gold data of supernova observations, it is 
demonstrated that this viscous cosmology model is rather well consistent to the observational data 
at the lower redshifts, and together with the diagnostic statefinder pair analysis it is concluded 
that the viscous cosmic models tend to the favored ACDM model in the later cosmic evolution, 
agreeable to lots of cosmological simulation results, especially to the fact of confidently observed 
current accelerating cosmic expansion. 

PACS numbers: 98.80.Cq, 98.80.-k 



I. INTRODUCTION 

Observations of type la supernova(SNe la) suggest that 
the expansion of the universe at later stage is in an ac- 
celerating phase. Additionally, the measurement of the 
cosmic microwave background (CMB) [jj and the galaxy 
power spectrum |2J indicate that in spatially flat isotropic 
universe, about two-thirds of the critical energy density 
seems to be stored in a dark energy component (the sim- 
plest candidate is the famous cosmological constant A) 
with negative enough pressure Q. Ironically, we do not 
know much about dark energy (DE) properties, if not less 
than those on the mysterious dark side of the universe Q . 

In order to explore the implying accelerating mecha- 
nism, many authors propose a variety of models to de- 
scribe the evolution of our universe, like the modified 
gravity [j| for example. Among these and opposite to 
the extending Hilbcrt-Einstcin action for general relativ- 
ity modifications, there exist a class of models that are 
based on searching for a proper equation of state(EoS) 
for the matter-energy fluid. Initially, this class of mod- 
els are exploited in the context of the perfect fluid. The 
viscosity concept is introduced into dark energy study 
relatively lately. And now it seems to play a more and 
more important and practical role in the more realistic 
cosmology model constructions. Other earlier attempts 
in this line can be found in references [a. YA UJl 12a . 133 . l33| . 
Additionally, for more details Gr0n have given a very use- 
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ful review for the subject in reference [Sj. 

Viscosity is a concept in fluid mechanics related to ve- 
locity gradient and is divided into two classes, shear vis- 
cosity and bulk viscosity. In viscous cosmology, shear 
viscosity comes into play in connection with space- 
time anisotropy. An analytic formula for the trace- 
less part of the anisotropy stress tensor has been de- 
rived by S.Weinberg^3|- Meanwhile, a bulk viscos- 
ity usually functions in an isotropic universe. Under 
the Friedmann-Robertson- Walker (FRW) framework, the 
energy-momentum tensor at most has a bulk viscosity 
term as = (p — (O)gik (( is bulk viscosity, and is the 
expansion scalar). Additionally, bulk viscosity related to 
the grand-unified-theory phase transition |ll| may lead 
to explain the cosmic acceleration expansion. 

At present, a large number of models exist but with- 
out very effective methods either verifying them or ruling 
them out. For this reason, there is a strong need for diag- 
nostic techniques. The statefinder diagnostic pair {r, s}, 
which is purely geometric quantities introduced by Sahni, 
Saini, Starobinsky and Alam[T3|. provides us a very use- 
ful method to discriminate cosmological models. The 
statefinder pair probes the expansion dynamics of the 
universe through higher derivatives of the expansion fac- 
tor 'a (a is the scale factor) which is a natural companion 
to the deceleration parameter that depends upon a. Dif- 
ferent models on the s—r plane accordingly show different 
trajectories. For example, the spatially flat ACDM sce- 
nario corresponds to a fixed point {0, 1} in the statefinder 
diagnostic pair {r, s} plane, with which the 'distance' of 
other Dark Energy (DE) models from ACDM can there- 
fore be established on the s — r plane [Tsj . Additionally, 
the statefinder pair has possessed the merits that can 
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discriminate among a large amount of models including 
ACDM, quintessence, kinessence, Chaplygin gas(see ref- 
erences El E Q E3) With the introduction of new 
observational techniques and increasing improvement of 
measurement, it is certainly for us to get more precision 
data and richer information about the geometric quan- 
tities, with possibilities to evaluate more practical cos- 
mology models, like those reasonable tries by considering 
viscosity to abandon the commonly used and simplest 
perfect fluid approximation to real cosmic media, and at 
that time some models would be either verified or ruled 
out by the statefinder diagnostic pair and other astro- 
physics observations. More reference about recent and 
future experiments can be found in the review paper [lr3|. 

This paper is arranged as follows: In Sec[n] the gen- 
eral formalism is presented for following discussions; In 
Scc lIIII we discuss the difference between perfect fluid 
models and non-perfect fluid models from the viewpoint 
of density and trajectories of the diagnostic pair{r, s}, 
especially in Quiessence model@,0,^3i I n Sec lIVI our 
attentions are focused on variable bulk viscosity dark en- 
ergy model with general EoS. The entropy of system is 
deeply discussed there apart from statefinder pair diag- 
nostics, as we have realized that the thermodynamics can 
reflect more important global characters for a compli- 
cated system, like our observable universe as shown by 
Brevik, Nojiri, Odintsov and Vanzoe in the reference |32j . 
And the last section is devoted to our conclusions. 



II. GENERAL FORMALISM 

Now, we introduce the basic framework for our discus- 
sions. That is, in Friedmann-Robertson- Walker(FRW) 
cosmology the metric of the system is chosen as: 



ds 2 = -dt 2 + a 2 ( 



dr 2 



1 - kr 2 



+ r 2 dVL 2 



(1) 



where a, and k are the scale factor and space curvature, 
respectively. The Einstein equations take the usual form 

G>„ = Rfj, v — -Rg^u = 8ttGT^ (2) 

Note that we have included the cosmological constant A 
in the energy-momentum tensor T^ v . In the FRW cos- 
mology with bulk viscosity, the stress-energy-momentum 
tensor is written as: 



{p + p)U^U v + pg^ - (911^ 



(3) 



where £ is the bulk viscosity, 8 the expansion factor de- 
fined by — Uft = 3ct/a, and the projection tensor 
is defined by h^ u — Q^y + UfJJ v with being the four 
velocity of the fluid on the comoving coordinates. On the 
thermodynamical grounds, £ is conventionally chosen to 
be a positive quantity and may depend on the cosmic 
time t, or the scale factor a, or the energy density p, 



etc. Through a series of calculations, the non-vanishing 
equations in Eq.(J2J are 
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P 



-3H(p + p) 



(4) 
(5) 



where p is an equivalent pressure defined by p = p — C,9 
and H denotes the Hubble parameter. Additionally, we 
use the unity convention 87rG/3 = c = 1. Normally, if 
we have known the information about both the equation 
of state(EoS) and the bulk viscosity £, the fate of the 
universe would have been determined by the Eqs.Q and 
©■ 

In the following parts, the viscous cosmology will be 
checked with the use of statefinder diagnostic pair {r, s} 
parameters. Here we first present their general expres- 
sions explicitly. The statefinder pair {r, s} is defined 
by (see reference |l2|) 
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where q = —d/aH 2 is the deceleration parameter. Com- 
bined with Eqs.Q and (0, Eq.ljlJ) becomes: 
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From the above formulas, we can see that the diagnostic 
statefinder pair r and s are related to the quantities p, 
p, p, £ and among which p, p, p can be reduced to one 
quantity if we have known the EoS. And for the little 
known global quantities £ and £ in the viscous cosmol- 
ogy, we can discuss the models with the simplest case 
of C = Co (constant) first as shown in the next section 
Recently, some authors have proposed a few fresh opin- 
ions about the possible forms of bulk viscosity £ such 
as £ = t9 for discussing dark energy cosmology and dark 
fluid properties (see reference ^3,E3); which has builded 
up a relationship between the bulk viscosity £ and the 
scale factor a. We will discuss such models with variable 
bulk viscosity in the section Hvl 



III. COSMOLOGY WITH CONSTANT BULK 
VISCOSITY 

In this section, we treat the universe model with a more 
realistic situation as containing two main media parts , 
that is, one is the mainly non-relative matter compo- 
nent M while another is the dominated dark energy X 
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component. Thus, the total pressure and density can be 
expressed as 



V = PM + PX — Px, P = Pm + Px 



(8) 



where the pressure of matter pm is a negligible quantity. 
And the EoS of the X part can take a usual factorization 
form px — wpx{w is called state parameter). 

From the viewpoint of the bulk viscosity, the simplest 
case is thought to be a constant bulk viscosity C = Co- In 
this section we mainly discuss the statefinder diagnostic 
pairs to two well known cosmological dark energy mod- 
els, added with a constant bulk viscosity. Through the 
comparisons of viscous and non- viscous models, it is ben- 
eficial for us to understand the role of cosmic viscosity, 
the properties of the cosmic models with common EoS 
and further our physical universe more comprehensively 



A. ACDM model 

So far as we know, the ACDM model, with mainly two 
components: the cosmological constant and cold dark 
matter, may be the simplest and the most consistent one 
with observation data. And also it has been deeply stud- 
ied with no viscosity assumption. In this subsection, we 
will continue to study it in the context of the viscous 
cosmology by using the statefinder diagnostic pair. 

Considering the Einstein equation |0J and energy con- 
servation equation (|3J), the following integration is obvi- 
ous: 



dp M 



Pm - 3CoVPM+PA 



da 



(9) 



where px has become p\ (the vacuum energy density). 
By working out the above integration, we get the relation 
between p and a: 



(P - Ko\/p - Pa) 



B 



where the B is an integration constant and the £ is a 
constant defined by 



3Co 



V^PA+MI 



Figure n demonstrates the a(t)) — p relation with the dif- 
ferent values of Co- There is a minimum value of the p 
denoted by p e j[ called the effective vacuum energy den- 
sity(EVED) corresponding to a — > oo . The expression 
of EVED is as 

p^i^Co + vV + gCo 2 ) - (10) 



The equation (|10|l has such a limit clearly: 

Pa 



lim p e ^ — 




FIG. 1: On the a — p plane, we can see the differences at 
the different values of £o- The horizontal lines represent the 
values, to which the next higher curves approach. Here we 
assume for simplicity the parameters(pA = 1 and B = 1). 
Note that £o = means without viscosity. 



which returns to the non-viscosity situation directly. 
Since pa denotes the vacuum energy density under the 
non-viscosity case, we might as well call it the conven- 
tional vacuum energy dcnsity(CVED) in contrast to the 
EVED. 

For the simplicity to discuss the trajectories on the 
figure which corresponds to the different values of the 
Co, we apply the relative changing rate as ^ defined by 



Ap 



p(a) 



\p{a + Aa) - p(a)\ 
p(a) 



(11) 



where the A a denotes a small change of the scale fac- 
tor a(t). After the interference of the bulk viscosity, the 
densities as a whole increase much more while the Ap 
changes a little from numerical analysis and the figure ^ 
and then the <j becomes smaller than before. And the 
trend of the changes is: the larger the bulk viscosity (Co), 
the bigger the value of EVED and the smaller the <;. We 
may sentence that the bulk viscosity stabilizes the evo- 
lution of the density and blocks then the rapid change of 
the universe, which is easily understandable by physics 
intuition as in the friction case that the friction force 
hinders dynamically the rapidly kinematic movements. 

Under this model, we put p = pu + PA, P = Pm + PA 
into Eq.(|7J) and get the statefinder pair as: 
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(12) 
(H>) 

do 



where the is the vacuum density parameter defined 
by r» A = pa/H 2 . 
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FIG. 2: On the s — r plane, we can see the differences with 
the different values of fo • When (o takes a smaller value, the 
trajectory looks not similar to a parabola. However, as £o be- 
comes bigger and bigger, its trajectory tends like the parabola 
more and more. Here we assume the parameter(pA = 1) for 
simplicity. Note that £o = 0(a circle) means no viscosity in 
this case. 



The expressions (|12|) demonstrate that the diagnostic 
pair {r, s} is dependent on a single dynamically changing 
quantity H. We draw on the figure [3 the trajectories of 
the statefinder pair {r, s}. The point {0, 1} corresponds 
to the diagnostic pair {r, s} for the ACDM model with 
no viscosity, and the curves which with the larger bulk 
viscosity £o are more parabola-like represent viscous sit- 
uations. The ACDM models that either are viscous or 
non-viscous are therefore differentiated by the trajecto- 
ries on the s — r parameters plane. 



B. Quiessence model 

The Quiessence model which characterizes itself with 
the EoS : 

Px = w p x 

(wq is constant, but asked not -1 as reason shown below, 
with contrast to the cosmological constant case) has been 
used recently to describe the dark energy behaviors. It 
is interesting to consider such model under the viscous 
situation. After assuming no interaction between two 
fluids M and X, we can decompose them , that is, the 
X and M satisfy the Eqs.Q and JSJ), respectively. By 
writing the X part out independently, we have: 



dp 



x 



(1 + w )pj 



= -3 



^ (13) 
a 



By working the above integration out, we get 

2 

Jt,0 - 3 ' 1 + m o) 

PX 



-^ + B(l + z)^ 
1 + w 



FIG. 3: In the p — z plane, we can see the differences at 
the different values of £o . The viscosity makes the density 
bigger than that of no-viscosity cases. Here we assume the 
parameter(_B =(7 = 1, and wo = — 0.5or — 1.5). Note that 
Co = means that there is no viscosity. 



where the redshift is denoted by z with z — ao/a — 1 
and B is an integration constant. And the density of 
the non-relative matter M component has possessed the 
following scaling relation from Eqs.Q and JSJ) 

p M = C(l + zf 

with a positive integration constant C . Then the value 
of the total energy density as an additive quantity can 
be easily composed as: 



3C ° +B(l + z)^ 



1 + w 



+ C(l + z) 3 (15) 



(14) 



It is a function about the scaling relation with variable 
z. Moreover, the concept of relative changing rate pa- 
rameter $ will be used again to investigate the effects of 
the bulk viscosity for clarifications. 

For the simplicity of discussing the parameter we 
first draw the relation of the Eq. i|15fl on the figure 
The definition of the <; is modified in this case as 

= \p(z + Az)-p(z)\ 

where the Az denotes a small change of the redshift z. 
With the same analysis as in the above subsection IIII Al 
we can still obtain the conclusion that the bulk viscos- 
ity stabilizes the evolution of the density in this model. 
It is worth noting that on the figure the trajectories 
are divided into two classes corresponding to wo > — 1 
and wo < — 1 respectively, and wo — — 1 is correspon- 
dence to the singularity as shown in Eq.(15) clearly. The 
phantom dividing phenomenon|l9j| can also appear on 
the statefinder pair plane as illustrated in the following 
discussions. 
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TABLE I: Relationship between Figures and parameters 





wo > - 


1 


w = —1 


Wo < — 1 


Co > o 


curves on 


4.2 6 


bullets on 4.1°,4.2 6 


curves on 4.1" 


Co = 


curves on 


4.3 C 


bullets on 4.3 C 


curves on 4.3 C 



bullet hereafter refers to point for diagnostic pair 
(s,r) = (0, 1), i.e., corresponding to ACDM model 
"4.1 denotes the first panel in the figure^] 
6 4.2 refers to the second panel in the figure 0] 
c 4.3 represents the third panel in the figure 0] 



The statefmder diagnostic pair of the Quiessence model 
can be gotten, when we put p = pu + Px and p ~ wopx 
into Eq.(0. The results are: 



3 o 9 Co , 
9 

r = -(1 + w )w £l x 



27 Co 
4 H 



(1 + 2w )w Q n x _ 3 Co 
2(w n x -3%) 2H 



1 

f 2 



x 



H 



(16) 
1) + 1 

G3>) 



where the fix denotes the dark energy density parameter 
denned by Q x = Px/H 2 . Considering Eqs.(QJ and (JT^J, 
we can ultimately transform q, r, s into such quantities 
depending on the rcdshift z only. 

The situation of the wo = — 1 is thought as the border- 
case called 'phantom divide' (see reference |21|). The 
singular property of the 'phantom divide' can also be de- 
scribed by the trajectories on the statefmder pair plane 
in the figure 0] The corresponding relationship between 
the trajectories on the figure0]and the values of the pa- 
rameters are arranged in the table [I] 

Obviously the point of (0, 1) that denotes the situation 
of the wo — — 1 divides the trajectories into two parts, 
the curves of Fig. 4.1 and the curves of Fig. 4. 3. The di- 
rections of the evolutions demonstrate that our universe 
is approaching the state of the wq = — 1 and then keeps 
itself stable there, which is consistent with results from 
lots of data analysis and cosmic simulations as favored 
to the ACDM model for describing later stage universe 
evolutions. 



in Ref. _20]), are made first, and then the trajectories of 
the model on the statefmder diagnostic pair plane are 
shown. For simplicity, we only discuss the universe filled 
by the fluid of only one component with the following 
EoS 221 



P= (7- l )p + Pi 



(17) 



where the p\ is a constant (as subscript often indicating 
present value). For system only consisting of same mass 
particles, like baryons, its pressure is proportional to nm 
where n is the number density with m as the mass for one 
particle. The 7 is the adiabatic exponent or barotropic 
factor, specifically, 7 = 5/3 in the case of extreme non- 
relativity, and 7 = 4/3 in the case of extreme relativity. 
As for other general cases, 7 will take complicated forms. 

The 7 — 1 will be equivalent to the state parameter w 
(hereafter) if the absolute value of p\ is much smaller than 
that of pressure p. The bound on the state parameter w 
is given out as (see reference [13,13, EH) 

-1.38 < w < -0.8 

and thus the adiabatic exponent 7 gets its own constraint 
as -0.38 < 7 < 0.2. 

By directly calculations from the Eqs.QJ and (JSJ, the 
scale factor is obtained as: 

Provided 7 < 0, 



a(t) = a 
If 7> 0, 
a(t) = a 



cosh 



2Ti 



7#o Ti sinh 



t-k 

2Ti 



cos 



t - tc 



O01T2 sin 



t - tc 



2T 2 J u ' V 2T 2 

where Ti,and T2 are defined by 

1 _ 1 



(18) 



(19) 



T 2 = 



And where 



7 = 7 - 3Ci 



IV. DE MODEL WITH VARIABLE BULK 
VISCOSITY AND EOS, THE ENTROPY 

In this section, we mainly discuss the viscous dark en- 
ergy model which characterizes itself with the variable 
state parameter w. We here take the form of the bulk vis- 
cosity as C = Ci<V a ( see reference 0,^|) which is a par- 
ticular situation of the general approaches as proposed 
in the reference |29|. Some thermodynamic discussions, 
especially the entropy expression (an instructive way to 
obtain the entropy expression for viscous cosmology from 
generalized Cardy-Verlinde entropy formula can be found 



is an effective parameter for the 7, and the po, 9q, to are 
the present energy density, expansion scalar, and cosmic 
time, respectively. 

Eas. (|18|) and i|19|) are compared by the use of curves on 
the figure |3 As for the case of 7 < 0, it can successfully 
explain the cosmic accelerating expansion in the late evo- 
lution universe. Conversely, it is impossible to do so for 
the case of 7 > 0, partly because there are negative val- 
ues of the scale factor a which seem un-physical by the 
direct mathematical treatment. So it is proper for us to 
only consider the case 7 < in the following discussion. 

The number of effective parameters from (|18H is re- 
duced to three: Ti, 9q, and 7. For the convenience of 




FIG. 4: The above three panels are numbered from up to down as 4.1,4.2 and 4.3: 4.1 is at the situation of wo = —1.5; while 
4.2 is in the case of wo = —0.5 and 4.3 is for viscosity free £o = condition. Arrows represent the directions of the evolutions 
of statefinder diagnostic pair about time. Here we have assumed for simplicity the parameters (B — C — 1) in the Eq. ljT^J . 



following discussions, it is beneficial here to consider pos- 
sible constraint to these parameters from the favorable 
fact a > particularly in the late universe evolution. By 
calculating a with the use of Eq.^HJ, we can have 



where W is defined by W = 76*0 T\ and x is defined by 
x = ao/a = 1 + z[z denotes the redshift). We have 
known that a > 0, 7 < 0, and T\ > 0, and thus a > is 
equivalent to 



37T1 



,37/2 



sinh(- 



2Ti 



'o^ , ,„ , ,i — to 
'■Wcosh(— — ) 



2Ti 



(20) 



W < 



1 - cxp(^) 
l + cxp(^) 



(21) 
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FIG. 5: On the t-a plane, the scale factor a(t) at the case 
of 7 > shows an oscillating expansion cosmology model 
while the 7 < case corresponds to the late time accelerating 
expansion universe 



Corresponding to t — > 00, as our interested later stage 
evolution for the observable universe, the right hand of 
the above inequality has possessed a limit as following 



-1 < 



1 - cxp(^fa) 
l + cxp(^) 



So the constraint l|21|) of W becomes 
W < -1 

Here we also give out, respectively, the expressions of 
density and expansion factor depending on the relevant 
variables t and z, which will be used in our entropy ex- 
pression calculations: 



9{t) 



po 



(1 



Pi 



Pi 



0(1 + *)^ 

Pol Pol 



(22) 



sinh(yf) + Wcosh(y^) 



W cosh( 



2Ti 

i-tn 
2Ti 



M^sinh(^) 



) = -4 ■ \Jl + (W 2 -l)(l + z)*r (23) 



(z) = ~w 

For the perfect fluid models in a closed cosmic system, 
the cosmic media is regarded without dissipation and the 
entropy is a conservation quantity with dS/dt — 0(5 de- 
notes the entropy of the system per unit volume) . How- 
ever, considering non-perfect fluid models, the entropy 
will change. Now we turn our attentions on the entropy 
of the model as introduced in this section. 

The relevant general formulas to be employed (see ref- 
erence UlSii) are: 



- T CT 



^e 2 -—( ),,(?• (2-D 



1 

kT 2 



where the is the entropy four-vector, r\ the shear vis- 
cosity, T the temperature, £ the bulk viscosity, <t^„ the 



shear tensor, 9 the expansion factor, k the thermal con- 
ductivity and as the space-like heat flux density four- 
vector. 

The entropy four- vector is defined by 



(25) 



where the na is the ordinary entropy per unit volume (n 
denoted the particle number per unit volume with a as 
the entropy of one particle). The expansion tensor 9^ is 
defined as: 

= -^{U Ka hZ + U v - ta h") 



The scalar expansion factor is 
tensor is defined as 



0{f = Uff,. The shear 



which is traceless, that is — and where the 
has been defined by h^ v = + V^U V . Defining the 
four-acceleration of the fluid as — — U^U^, the 
space-like heat flux density four- vector is given by 



Q" 



-Kh» v {T v + TA v ) 



In the case of thermal equilibrium, Q M = 0. 

Under the background of FRW metric, we can have 



fujy = 0, 



= 3H 



Eqs.|(23 and © yield 

S° = n<j, 5* = 0(i = 1, 2, 3. 



After taking account of FRW metric, we can get the fol- 
lowing differential equation as 



(26) 



where the ',0' denotes time derivative d/dt. We assume 
that the fluctuation of temperature is so small that it is 
negligible. The Ea. i|2fi|l can be transformed as 



^e 3 

3T 



(27) 



Then, based on Ea. (|18(l . we work the differential 
Ea. (|2*fj|) out and get: 



S° = 



l-W 2 9 2 /9 2 \ 69 Q T 



+ Ce w 



(28) 



where the C is an integral constant. From the equation 
of J2BJ, we find that when 9 = 9 /W, the S a approaches 
to infinity if we did not consider the values of the con- 
stant C . To avoid this un-physically mathematical phe- 
nomenon, we choose the value of C as 



C 



6 W 

6TW 3 
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which will lead to such a limit that 



lim S° 



3TW 



without the formal singularity. The expression l|28ll of 
the S° therefore becomes: 



S°(0) = 



' 6w 3 t ' w 2 e 2 - el 



or in a more easily reading form 



S°(z) 



2 + (W 2 - l)(l + z) 



37 



(29) 



(30) 



6W 3 T 

The conservation equation for the particle number is: 

(nLP) ; „ = 

which means that na 3 = constant in the comoving frame. 
Therefore the entropy for the whole observable universe 
in this model is 



S(z) = (na)a 3 = (^fS°(z) 
By considering Eas. (|23|) and l|28|l . we can obtain 



(31) 



S(z) 



~ 6W 3 T ' \l + z 
2 + {W 2 - + z) 3 ^ 



(32) 



With the above expressions, we have established the 
relation between the entropy S and the redshift z. 

Now, we will use 157 gold data points as presented by 
Riess et al|3lj to confront with our model and constrain 
parameters W as well as 7. In order to maximize the 
following likelihood function (see references p9l l30| ) : 



L cx exp[ — - ] 
we minimize \ 2 which here is expressed as 



X 



E 



(T, 



(33) 



(34) 



where 2^, ^th(z.i), and Oi are known from Gold data. 
jJLth{z) is defined by 



Hth{z) = 51g 



(1 + *) 



dz' 
hW) 
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where h(z) is the reduced Hubble parameter with h(z) = 
H(z)/H = 9(z)/8 n . Here we have assumed curvature 
k = and c is a constant ( also one of constrained quan- 
tities). 

Through the numerical calculations, we find that the 
best consistent values can be taken as 



W 



-1. 



7<0 or W < — 1, 7 = 




-0,2S -0.2 



0.15 -0.1 -COS 



FIG. 6: The likelihood contours obtained for W and 7 with 
157 gold SNe la data analysis. 



45 



i 40 



35 




+ gold data 
theoretical curve 



0.5 1 

z 



1.5 



FIG. 7: The comparison of experiment data and theoretical 
curve by fitting the 157 gold SNe la data 



Because the 7 has also possessed the constraint of 7 < 0, 
we can merely take W = -1 then. The 68.3%, 95.4% and 
99.7% likelihood contours are shown on the figureEl The 
next figure is the comparison between experiment data 
and the theoretic model estimations, from which we can 
see that the theoretic estimations can well fit the Gold 
data for the smaller redshifts. 

Taking W = — 1 into Ea. (|32fl . the expression of entropy 
is reduced into 



S(z) 



3T 



l + Z 



(35) 



in which the entropy density S° does not alternate, but 
the entropy S changes as the "volume" a 3 varies. Ob- 
viously the entropy of the Ea. (|35|l provides an arrow of 
time for cosmic evolution with the meaning that the en- 
tropy of our observable universe is always increasing. 
For the definition of the parameter W, W = ^OqTi , 
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TABLE II: Relationship between models and figures 



FIG. 8: On the s — r plane, the trajectory is a hyperbola. 
The present point is at r = l,s = 0. The arrows indicate 
the cosmic parameters changing tendency, which prefers the 
ACDM model in the later cosmic evolutions. 



parameter 7 when one takes W = — 1 becomes 

9pi 



El 

Po 



7 = 7 ^ 3& 
Taking the result into Eq.J23), we can get 



P = Pa 

which is a constant. Consequently, pressure p also takes 
a constant value. Comparing these results with the well- 
known ACDM model, we may conclude that in the vis- 
cous cosmology cases, the best fitting results still favor 
the ACDM model. This point can also be reached from 
the view point of of statefinder diagnostic pair. 

Via using the statefinder diagnostic pair {r, s} expres- 
sions, we can obtain directly, 



1^ 

§V(r- 



1+1 



^(7-1) 
V-l 



(36) 

jut) 



where V is defined by 



V 



Pi{ 

P Pa 



The trajectories is show on the figure |H1 When taking 
the best fitting parameter W — —1, we obtain V = 0, 
r = 1 and s = that is also the case of the ACDM 
model. So it confirms our previous point again that in 
the viscous cases, the best fitting model still returns to 
the conventional ACDM model. 



the s — r plane 


C = Co 


C = (ia/a 


ACDM 


FigH 


®° 


Quiessence 


FigH 


0" 


Variable DE EoS 
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Figun 



°Cg> denotes no figures 



V. CONCLUSIONS AND DISCUSSIONS 

We mainly discuss the behaviors of some viscous cos- 
mology models on the statefinder pair plane on the pur- 
pose to mimic dark energy characters, with the hope to 
demonstrate that cosmic viscosity can also play the role 
as a possible candidate for dark energy. To this aim, we 
are first to give out the formulas for viscous cosmology 
statefinder pair expressions J7J. After introducing the 
relative changing rate ?, the global evolution of density 
is found to relax as more stable in the viscous model sit- 
uations. At the same time, the trajectories of viscous 
universe statefinder pair on s — r plane become quite dif- 
ferent from non- viscous cases as tabic ILT1 

Particularly for the Quiessence model, the singular 
property of the 'phantom divide' w = — 1 can be clearly 
demonstrated on the statefinder diagnostic pair plane 
by the completely different trajectories to discriminate 
themselves. And the directions of the evolution of the 
trajectories on the statefinder diagnostic pair planes for 
the three models all point to the point of r = l,s = 
(ACDM model preferred), that is the universe favors 
simple in the later evolution stage with scale largely ex- 
panded. The ACDM cosmology is simply consistent with 
the current astrophysics observations, especially the cos- 
mic late time accelerating expansion, but the cosmologi- 
cal constant has been puzzling ever since. 

Additionally, we also have a try to describe the en- 
tropy of the viscous cosmology system. After adopting 
the EoS lfT7|) and the bulk viscosity £ = (iH, we deduce 
out the concrete expression for the entropy. Then the 
157 gold data from the supernova observation are used 
to constrain parameters W and 7, and we therefore get 
the most favorable parameter: W = — 1. Further we find 
that the entropy of the universe is always increasing with 
cosmic evolution, which is consistent with the thermody- 
namics arrow of time. 

Observational cosmology across this century has chal- 
lenged our naive physics models, and with the anticipated 
advent of more precious data we have the chance to un- 
derstand or uncover the universe mysteries by more prac- 
tical modelling. Quite possibly we will get more hints to 
unveil the cloudy cosmological constant puzzle. In the 
simple constant bulk viscosity case (a proto type or a toy 
cosmic media model) as demonstrated in section three 
the vacuum energy density can be shifted by the bulk 
media viscosity to arrive at an effective vacuum energy 
density (EVED) or we may say that the constant viscos- 
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ity can tune the cosmological constant in a sense if we 
have possessed a suitable cosmic media model. We ex- 
pect more encouraging work on non-perfect fluid cosmic 
concord models to come soon and we believe this line of 
trying can contribute us new understandings to the mys- 
terious dark side of our complicated but observable and 
conceivable universe. 
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